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We study dynamical heterogeneity and glassy dynamics in a kinetically constrained lattice gas 
model which has both translational and rotational degrees of freedom. We find that the rotational 
diffusion constant tracks the structural relaxation time as density is increased whereas the trans- 
lational diffusion constant exhibits a strong decoupling. We investigate distributions of exchange 
and persistence times for both the rotational and translational degrees of freedom and compare our 
results on the distributions of rotational exchange times to recent single molecule studies. 

I. INTRODUCTION 

Glassformers are dynamically heterogeneous. Neighboring regions, nanometers in size, can have local relaxation 
times which differ by several orders of magnitude 0, Q • Experimental measures of heterogeneous dynamics often 
probe rotational degrees of freedom. For example, Deschenes and Vandenbout [3, |3j have measured the rotation of 
single probe molecules in a polymer film near the glass transition. Here, we investigate a kinetically constrained 
lattice gas model with translational and rotational degrees of freedom which captures many of the essential features 
seen in these experiments such as stretched exponential decay of rotational autocorrelation functions as temperature 
is decreased and heterogeneous distributions of exchange times. 

In section HT1 we present the model and the computational methods used. Section ITTT1 demonstrates the existence of 
heterogeneous dynamics. We also present ensemble measurements which display a precipitous dynamical slowdown 
and breakdown of mean field dynamical relations. Section llVI measures distributions of persistence and exchange 
times. Finally, section compares these distributions to those observed recently in single molecule studies. 

II. MODELS AND COMPUTATIONAL DETAILS 

One route to understanding dynamical heterogeneity relies on the presence of local steric constraints on the move- 
ment of particles which make themselves felt to an increasing degree as temperature is lowered (or density is increased) . 
The kinetically constrained lattice gas models Q are simple caricatures of glassformers which employ local steric con- 
straints as their sole means to glassiness in the absence of any non-trivial static correlations between particles. It 
has been shown that, despite their apparent simplicity, these models exhibit surprisingly many of the hallmarks of 
glassy behavior that have been the focus of recent experimental and theoretical efforts [J 0, |(J 0, @- The model we 
study consists of hard core particles on a triangular lattice with no static interactions other than those that prohibit 
multiple occupancy of a single site. To each particle is associated a vector which can point along the six directions 
bisecting the triangular lattice. In other words, the vectors point toward the interstitial sites of the lattice 0. 

Translation of particles obeys the kinetic constraints of the two vacancy assisted triangular lattice gas, or the (2)- 
TLG a particle at site r is allowed to move to a nearest neighbor site, r', if (1) r' is not occupied and (2) 

the two mutual nearest neighbor sites of r and r' are also empty. These rules coincide with a physical interpretation 
of steric constraints on the movement of hard core particles in a dense fluid |Toj . 

Rotation of particles obeys a similarly physically motivated kinetic constraint if one imagines that the particles 
have small hard protrusions along the direction of their orientation vector 0: a particle, i, with a unit vector p^, can 
rotate either 60 degrees clockwise or counterclockwise provided the two neighboring lattice sites along the direction 
between the initial and final orientations of Pi are empty. Translations preserve the direction of p^. We refer to this 
model as the rotational TLG 9] . Due to the absence of non- trivial static interactions, there are no static correlations 
between particles. However, constraints on the kinetics allow non-trivial dynamic correlations to emerge in trajectory 
space. 

In the computer simulations, we investigated particle densities, p, between 0.10 and 0.81 on a lattice with edge 
length L = 128 (L = 256 for p = 0.81). The density p — 1 corresponds to the completely full lattice. At densities 
up to and including 0.77, over 60 independent trajectories of lengths 10-100 times r a , where r Q is the time for the 
self-intermediate scattering function at q = (7r,0) 01 to reach 1/e of its initial value (see below), were run. At 
densities 0.79, 0.80 and 0.81, four to sixteen trajectories were run. These trajectories were stored logarithmically for 
later analysis (i.e. configurations were saved after 1, 2, 4, 8, 16, 32, etc. sweeps). Time was measured in Monte Carlo 
sweeps. During each sweep, particles were chosen randomly and translational and orientational moves were attempted 
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FIG. 1: Rotational trajectories of single particles at (a) low density (p = 0.30) and (b) high density (p = 0.77). 



with equal probability. For the higher density runs, a continuous time algorithm was used for greater efficiency 13] . 
This algorithm involved making and updating a list of only those particles which have the possibility of either moving 
or rotating and choosing from among those exclusively during every move. The time was then incremented as one 
over the number of possible moves. Due to the lack of static correlations, initial configurations could be generated by 
random occupation of empty lattice sites by particles with random orientations until the desired density was reached. 

III. HETEROGENEOUS DYNAMICS AND ROTATIONAL CORRELATION TIMES 

Despite the lack of static correlations between rotational and translational degrees of freedom, dynamic coupling 
exists via the kinetic constraints 0- A particle which cannot rotate because it is sterically blocked by neighboring 
particles must wait for those particles to translate away before it is allowed to rotate again. In this sense, we expect 
rotational dynamics to be a good indicator of local structural relaxation in this model, as they are in experiment. 

Fig- HI shows rotational trajectories of six particles at low and high densities. The particles at low density (Fig. QJi) 
rotate freely, performing a random walk through all angles. At higher densities, the particles perform random walks 
punctuated by periods of little to no rotation. In contrast to the trajectories at low density, dynamics at high density 
are clearly heterogeneous: at any given time, some particles are rotating quickly while others are essentially frozen. 
Similar behavior has been observed for translational motion of a probe molecule immersed in spin-facilitated models 
[Hill and in poilymcr films 0,0. 

A useful ensemble measure of slow dynamics is the rotational autocorrelation function, C r (t) = (Pi(0) ■ Pi(t)), 
which indicates the time it takes for a particle to lose memory of its initial spin orientation 2, 9J. Here, the angled 
brackets denote an average over all particles and times, t. A plot of C r (t) is shown in Fig. [2^. At low densities, 
relaxation shows a simple exponential profile. As density increases, the curves become more and more stretched 
exponential indicative of averaging over multiple relaxation timescales. This stretched exponential decay is what one 
would expect qualitatively from rotational trajectories such as those depicted in Fig. ^p. An important quantity that 
can be extracted from C r (t) is the rotational correlation time, r r , which is defined as C r (r r ) = 1/e. The inverse of 
this timescale is the rotational diffusion constant, D r — t~ . 




FIG. 2: (a) Rotational autocorrelation at (from left to right) p = 0.50, 0.60, 0.65, 0.70, 0.75, 0.77, 0.79 and 0.81. (b) Scaling 
of the rotational diffusion constant, D r (stars), and translational self-diffusion constant, D s (open circles), with structural 
relaxation time, r a . The dashed lines are the power laws, r" 1 and r~ ' 6 as indicated. 

Translational relaxation is often studied via the self-intermediate scattering function, F s (q,t) = (e lq '( ri (*) -ri (°))}. 
Here, rj(i) denotes the position of particle i at time t. The decay of the scattering function to 1/e at wavevector 
q = (?r, 0) is typically defined to be a structural relaxation time, t q , as it gives a sense of how density fluctuations 
relax at relatively short lengthscales. Dynamic behavior at large lengthscales is studied via the self-diffusion constant, 
D s , extracted from the mean-squared displacement, (\Ar i(t)\ 2 ) = — 1^(0)1 2 ). The self-diffusion coefficient, D Sl 

is defined as D s = limt_ >rm ([Ari('t')| 2 )/4t. We omit further discussion of these quantities as they have been presented 
at length elsewhere for this 0, OH and other models Q . 

An important ramification of heterogeneous dynamics is the breakdown of mean-field dynamical relations such 
as the much studied Stokes- Einstein (SE) relation 0, 0, |^. In a system with homogeneous dynamics such as 
a normal liquid, we expect relaxation behavior to be similar at all but the smallest lengthscales. In a glass, the 
presence of dynamic heterogeneity implies that such mean-field relations can be violated. In Fig. [2Jd, we plot the 
rotational diffusion constant, D r , and the self-diffusion constant, D s , versus the translational structural relaxation 
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FIG. 3: Distributions of translational (a) exchange and (b) persistence times at (from left to right) p — 0.50, 0.60, 0.65, 0.70 and 
0.75. The inset to (a) shows the exchange time distributions as a function of linear time, t, as opposed to logarithmic time to 
emphasize the emergence of broader long time tails as density is increased. The two distributions are related via P x> trans (log t) = 
tP x ,trans{t) In 10. The solid lines in (b) are the distributions of persistence times calculated from the distributions of exchange 
times via eqn.Qand the open circles are results of direct calculations. 



time, r a . Rotational diffusion tracks structural relaxation whereas self-diffusion does not. That is, D r ~ t^ 1 whereas 
D s ~ Tq' 6 - The same trend has been seen in experiment [T^ |. The scaling of D r with r Q can be rationalized 
qualitatively from the idea mentioned above that the ability to rotate is intimately tied to local structure. 



IV. DISTRIBUTIONS OF EXCHANGE AND PERSISTENCE TIMES 



Direct measures of het erog eneous dynamics in glassy systems are distribution functions of persistence, P p (t), and 
exchange, P x (t), times |2l ll4lll8l ll9j. Persistence times measure the first instance of a change in state given an initial 
configuration. Exchange times measure the duration of particular states. For example, the distribution of rotational 
persistence times, P p ro t(t), in this model is the distribution of times, given an initial configuration, when a particle 
changes its rotational state for the first time. The distribution of rotational exchange times, P x ,rot(t), is the length 
of time a particle remains in a particular rotational state. These distributions are multi-point functions because they 
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depend not only on two points in time, but on all intervening points as well. Fig. |3i and Fig. Up, show the distribution 
of exchange times for translational and rotational motion, respectively. Distributions of exchange times are important 
for understanding the origin of dynamical decoupling phenomenon (see Fig. 03). 

The distribution of exchange times for translations and rotations are very similar. There is very little change in 
short time structure as p increases. Evidence of broadly distributed dynamics occurs in the long time tails. The 
distributions of persistence times for translation and rotation are also qualitatively similar and display the same 
features noted in previous studies 0, ITsI ] . Differences in structure between rotational and translational persistence 
can be attributed to the differing kinetic constraints. A particle can have the ability to rotate while being constrained 
translationally and vice versa. 

Recently, it has been shown that the distributions of exchange and persistence times are related via the equation 

/OO 
P x (t')dt'. (1) 

The constant of proportionality is fixed by normalization. This should be a general relation independent of the model 
or observable studied. The solid lines in Fig. [3Jd and Fig. 0]d show the distribution of persistence times calculated via 
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FIG. 5: (a) A plot of the absolute rate of angular change, |A#/At|, for a rotational trajectory of a single particle at high density 
(p = 0.75). Three examples of exchange times are indicated by the double headed arrows, (b) Comparison of distributions of 
rotational exchange times obtained via data like those shown in panel (a) (open circles) and directly (solid line, see Fig. HJi, 
inset). The inset shows a comparison of the distribution of rotational persistence times. The values for the open circles in the 
inset are obtained by applying eqn.^to the the data in the main part of panel (b). 

equation (JJ. The open circles are the results of direct calculations. We see that equation JIJ is well-satisfied. 

V. COMPARISON WITH EXPERIMENTS 



Recent single molecule experiments determined distributions of rotational exchange times by following the dipole 
of embedded dye molecules Near the glass transition temperature, the absolute value of the rate of angular 

change, |A0/At|, showed abrupt changes between different dynamical environments. We can anticipate from Fig. 

that the same quantity in the rotational TLG will show similar behavior. From this quantity, Deschenes and 
Vandenbout extracted a distribution of exchange times using a standard deviation criterion 0]. Here, we examine 
to what extent the distribution of times measured in this way corresponds to the distribution of rotational exchange 
times defined in section Hvl 

In Fig. G3i, we plot |A0/At| where AO = 9(t + At) - 6(t) for the rotational TLG at high density. Deschenes 
and Vandenbout assigned an exchange event whenever the average angle jump changed by more than two standard 
deviations from the previous average angle jump. Due to the coarse grained nature of the rotational TLG, \A6/ At\ 
changes by discrete jumps. Therefore, exchange events can be unambiguously assigned whenever such a jump occurs. 
Fig. |SJi plots \A6/At\ in increments often sweeps during a portion of a single molecule trajectory. The fastest exchange 
event measurable at this time resolution is 10 sweeps. With infinite time resolution, however, we see that the definition 
of exchange times used in section llVl would correspond precisely with the Deschenes- Vandenbout procedure. 

We verify this in Fig. [SJx Here, the distribution of times obtained via the Deschenes- Vandenbout procedure outlined 
above (open circles) [21j are compared to the distribution of rotational exchange times calculated in section llVl (solid 
line). The comparison is very good for about ten orders of magnitude. This method overestimates the distribution 
at early times for the reasons of time resolution mentioned in the previous paragraph. The inset shows a comparison 
of the distribution of rotational persistence times. The open circles are obtained from the open circles in the main 
panel of Fig. |3p via equation QJ. Once again, the agreement is very good. In particular, the data obtained via the 
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procedure of Deschenes and Vandenbout captures the structure and location of the main peak. 

Single molecule experiments also measured the mean rotational exchange times, T xrot (i.e. the first moment of 
Px.rot(t), inset), as a function of temperature. It was found that T x , rot scaled with temperature in the same 

way as r r , the rotational correlation time, and that T x ^ rot was 10-20 times larger than t> at a given temperature 
Our calculations indicate that T x ^ rot ~ (1— p)~ 2 or, using the definition that 1— p = exp(— 1/T) [||, Txj^nt ~ exp(— 2/T). 
The rotational correlation time, r r , scales like the structural relaxation time: lnr Q ~ cxp(1.7/T) [5J,la|- That is, t x rot 
exhibits a much weaker dependence on temperature than r r similar to what is found in spin- facilitated models |l9j |. 
We also find that T XiTOt is as much as six orders of magnitude times smaller than r r at the highest densities. This 
difference with experimental findings is most likely due to the issue of time resolution mentioned earlier and also 
discussed in [20] which biases experimental measures of r x rot to longer times. 
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